Finite Temperature Depinning of a Flux Line from a Nonuniform Columnar Defect 
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A flux line in a Type-II superconductor with a single nonuniform columnar defect is studied by a 
perturbative diagrammatic expansion around an annealed approximation. The system undergoes a 
finite temperature depinning transition for the (rather unphysical) on-the-average repulsive columnar 
defect, provided that the fluctuations along the axis are sufficiently large to cause some portions of 
the column to become attractive. The perturbative expansion is convergent throughout the weak 
pinning regime and becomes exact as the depinning transition is approached, providing an exact 
determination of the depinning temperature and the divergence of the localization length. 



The finite-temperature statistical behavior of magnetic 
flux lines (FLs) in Type-II superconductors with various 
kinds of disorder has been a topic of considerable inter- 
est and technological relevance following the advent of 
high-T c superconductors about a decade ago 0. Early 
on, it was realized that compared to point defects, ex- 
tended defects such as twin boundaries and columnar 
amorphous regions created by heavy ion irradiation |^] 
are much more effective in pinning of FLs at higher tem- 
peratures, an essential feature for technological applica- 
tions of these new materials. Most analyses of extended 
defects studied uniform structures However, it is nat- 
ural to expect that such defects contain spatial inhomo- 
geneities in themselves. This paper addresses the impli- 
cations of such a possibility for the simple case of a single 
FL in the presence of a single inhomogenous columnar de- 
fect with short-range correlated disorder along its axis, 
which is aligned with the external magnetic field. 

A number of closely related works on 2D wetting with 
quenched disorder Q^], directed polymers with random 
interaction j7j, strongly inhomogeneous surface growth 
||, and 2D wetting and directed polymers in a periodic 
potential |^|| have mostly considered a (1+1) dimen- 
sional geometry, for which the disorder was marginally 
relevant at the depinning transition. Higher dimensional 
systems were glossed over upon the observation that dis- 
order was irrelevant, i. e., the critical behavior near the 
depinning transition was the same as the uniform defect 
case. Nevertheless, it would be very useful to compute 
the change in the free energy and the localization length 
due to quenched disorder, in some vicinity of the depin- 
ning transition where small-scale details of the disorder 
becomes unimportant. This paper presents a controlled 
perturbation expansion, which becomes exact in the weak 
pinning limit, for the corrections to the free energy of the 
FL due to fluctuations in the pinning potential along the 
columnar defect. 

Consider a single FL (with line tension e) in the pres- 
ence of a columnar defect of radius 6, oriented in the 
z— direction with a pinning potential U (z) inside the de- 
fect. The Hamiltonian for a FL configuration {r(^)} is 



W({r(z)}) = 




[U Q + SU(z)} h(r(z)) 



(1) 



where h(r) = 0(6— |r|).(0 is the unit step function.) The 
fluctuations SU(z) around the mean potential Uo have a 
Gaussian distribution with zero mean and correlations 



5U(z)SU(z') = — exp 

27T 



{z-z'f 
2£ 2 



(2) 



where the overline represents averaging over the disor- 
der potential. Various thermodynamic properties of the 
system such as the free energy can be deduced from the 
partition function of the system: 



Z = Z: 



2?r(z)exp[-/3H({r(z)})], 



(3) 



where (3 — l/T (k^ = 1 throughout), and the normaliza- 
tion factor is chosen such that Z = 1 for Uq = 8U(z) = 0. 
Z is a random variable that depends on the particu- 
lar realization of disorder. Thermodynamic quantities 
typically involve averaging of In Z, which is usually too 
difficult to obtain directly. An alternative approach is 
to calculate moments of the random variable Z since 
Z n = e nhiZ is the characteristic function for the ran- 
dom variable InZ IllOtl. Indeed, 



Z n = exp 



n K 



£^ fe (lnZ) 



k=l 



(4) 



where Cfc(lnZ) is the fcth cumulant of \nZ. The nth 
moment of Z is 



Z n = J Pr 1 (z)...Pr„(z)exp|-/3 J dzj^ 
+ [Uo + SU(z)]J2h(r l (z)) \ . 



2\dz 



(5) 



Upon averaging over disorder, one obtains expectation 
values of the form 



1 



-0 J dzm(z)SU{z) _ Si J dzdz' rn(z)m(z')SU{z)SU(z') 



exp 



(?a\ 



dz [m(z)Y 



(6) 



where m(z) = X)"=i M r «( z )) ^ s the number of lines that 
are inside the defect at height z, and the approximation 
to(z') « m(z) leading to the final result is good provided 
that £ < £z = /3e6 2 , the diffusion length inside the defect. 
Using [m(z)] 2 = m(z) + m(z)[m(z) — 1], the nth moment 
can be rewritten as 

4 — 1 

/3 2 a 2 £ | dz^Mz))/^*)) . (7) 



x exp 

where 

Wo({r*(*)}) = 



I Iff = 



+ U eS h(r(z)) } (8) 



is the Hamiltonian corresponding to a FL in the presence 
of a uniform columnar defect, with defect energy per unit 
length 



U eS (T) = U - a 2 £/2T. 



(9) 



The solution to the uniform cylindrical defect prob- 
lem is well studied The free energy per unit length 
/o and probability density ip 2 (r) 01 the FL at a posi- 
tion r are given respectively by the ground state energy 
and wavefunction of the corresponding two-dimensional 
"Schrodinger Equation" 



r 2 



2c 



^Vi + U cS h(r) 



V'o(r) = Mo(r). 



For U e g < 0, the FL is pinned to the defect, with 



Jo 



l^effl (T\ 



-2(T/T t ) T > T 



dr r Vo( r ) 



1/2 



be {T/Te) \ T> T> 



(10) 

(11) 
(12) 



where £± is the localization length and Tg = \/ — U c ff£ b 
is the localization temperature, below which the FL is 
strongly pinned to the defect (£± ~ b) and the free energy 
is dominated by the ground state energy rather than the 
entropy (fo ~U cS ). 

For Uq < (on-average attractive defect), U s < 
and the FL is pinned to the defect at all temperatures. 
However, for Uo > (on-average repulsive defect), there 
is a depinning transition when U c s = 0, at a temperature 



T d 



2tV 



(13) 



For T > T d , the FL is delocalized, with f = 0. As 
t = (Td — T)/Td \ + , the free energy and localization 
length diverge as 



2 tat 



fo{t) « - — 7-4 exp 



2^ 
tat 



£±(t) « b exp 



la^ 

t at 



(14) 
(15) 



where a c = C/ 3/4 ? 1 / 4 (26/0 1/2 
Notice that 



exp 



/3 2 a 2 £ 



(16) 

where Zo = e _ ^ oi and (...)o denotes an average over a 
canonical ensemble of n non-interacting FLs in the pres- 
ence a uniform defect of strength U e g. This additional 
exponential factor is due to the fact that all n FLs sample 
the same "quenched" realization of disorder, and can be 
interpreted as an effective two-body interaction between 
FLs that is nonzero only when both FLs are inside the 
defect. 

For T > Tf,, the FLs are delocalized and the free energy 
for the quenched system is identically 0. In the weak pin- 
ning regime Td > T 3> Tg, corrections to the free energy 
can be computed as a controlled perturbation series in 
some "small" parameter. Expanding Eq.([l6|) in a power 
series: 

^ = 22 fc! 22-22 J 

x (ft( r<1 {z x ))h{T jl (z x ))...h{r ih (z k ))h(r jk (z k ))) Q . 

The expectation value factors for distinct FLs since each 
FL wanders independently of each other. After ordering 
in z such that z% < Zj for i < j, the product can be 
further factorized as 



2s TT (H r (zi))H r ( z i+i)))o 



mh(r(z t ))) =ku 

\i=l I »=1 

* 2 = (h(r(z))) = J dv^{v) 

|r|<6 



, (18) 



(19) 



since each visit to the defect (at z = Zj) is a renewal 
event, i. e., the configuration of the FL for z > z, is 
independent of its configuration for z < z.j. However, the 
probability of a subsequent visit to the defect is enhanced 
up to a distance of order l z = £, z (£±/b) 2 , after which the 
FL "forgets" its previous visit. This enhancement can be 
captured in an effective "propagator" 



2 



L — L Q exp 




(a) 



= X(z-z') 



X - 



T 2 



(b) 



FIG. 1. (a) Diagrammatic representation of the disor- 
der-averaged partition function, (b) Feynman rules for cal- 
culating various diagrams that contribute to vertex functions 
V k . 




(a) 





FIG. 2. (a) Diagrammatic representation of the geometric 
series that contributes to V%. (b) Leading-order contribution 
to V 3 . 



(h(r(0))h(r(z))) o - (fa(r(0)))g 
(h(r(0))} 2 



(20) 



which contains all the z— dependence and becomes expo- 
nentially small for z > £ z . Then, 

(f[h(r( Zi ))) =^/f[[l + X (z i+ i-z i )]. (21) 

\i=l / o i=l 

Upon further manipulation, the disorder-averaged parti- 
tion function can be recast into a form [See Fig. |l](a)] 



Z n = e - n foL/T exp 



n(n — 1) 



k=2 



(22) 



where the vertex functions T4 are a sum over all con- 
nected graphs with k incoming and k outgoing lines, con- 
structed from the building blocks shown in Fig. |l](b). The 
combinatorial factor coming from the choice of the first 
pair to be contracted is pulled out for notational conve- 
nience. 

Leading-order contributions to the free energy come 
from V2, which can be computed exactly in terms of \l/o 
and x- As shown in Fig. ||(a), the terms from a geometric 
series, and 



coefficients of terms in the exponent of Z n that are lin- 
ear and quadratic in n, respectively. 



^/"quenched 



2T d 1 - [C**K,/T*) 



5f? 



2Ca 4 . 



exp 



ta 4 



quenched 



^/quenched- 



(24) 
(25) 



Thus, the free energy is self averaging as L — > 00. Note 
that the geometric series converges when a 2 > \j2Ca 2 c . 
The leading-order contribution to V3 is [See Fig. ||(b)] 



(n 



dzx(z) 



(26) 



Near the transition, f™dzx(z) ~ (T/Ti) 4 /^ and 
therefore this next correction to ^ /quenched is smaller by 
a factor of order (a 4 /a 4 t) 4 exp(— 2a 4 /ta 4 ), in fact, lead- 
ing order contributions from higher order vertex func- 
tions form a series where each term is smaller from the 
previous one by this same factor. Thus, the asymptotic 
behavior of the free energy and the localization length 
near the depinning transition is given by Eqs.(|l4j) and 
(|l5|), since 



2 c 00 

i=0 



2 

-*n 



dz x 2 (z) 



(23) 



Just below the pinning temperature, < t <C 1, 
^0 Jo°° dzx 2 ( z ) ~ CCz j where C is a constant of 0(1). In 
this weak pinning limit, the correction to the disorder- 
averaged free energy and its variance are given by the 



lim 

t->t: 



quenched 

To 



0. 



(27) 



The perturbation series can be used to compute cor- 
rections to the free energy even away from the depinning 
transition, provided that r > Tj, i. e., within the weak 
pinning regime. The situation is clear for Uo < 0. In 
the interesting case Uq > and o o Cl this condition 
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FIG. 3. Top : Diagram depicting three possible scenarios 
as a function of a and Uq. Bottom : Depinning temperature 
Tb as a function of disorder, for a "repulsive" defect with 
Uq > 0. The diagrammatic series gives a good approximation 
to the pinning free energy whenever T>T. 

breaks down near a temperature T » (cr 2 /2cr 2 ) 2 ^ 3 T^. 
For T < T, the localization length becomes of order 6 
and energetic contributions dominate over entropic con- 
tributions to the free energy. At T — 0, the FL sits adja- 
cent to the defect and makes excursions into it whenever 
energy fluctuations are favorable. The energy (per unit 
length) cost of a typical excursion of length £ is 



E(eye « ib 2 /£ 2 + u - V^/i, 



(28) 



which is minimized for I* w (16e 2 & 4 /cr 2 £) 1//3 , and yields 
an estimate for the ground state energy 



f(T = 0) 



Uo 
2 



2a 2 



2/3 



(29) 



For a < er c , fluctuations are not favorable enough to make 
excursions into the defect, and the FL is depinned at 
all temperatures, which suggests why the perturbation 
expansion breaks down for this case. 

The T = behavior can be expected persist to finite 
temperatures as long as i* < £ z , or equivalently, T > 
T . Thus, as temperature is increased, there is a simple 
crossover from strong to weak pinning near T. Fig. |3| 



depicts a cartoon phase diagram and the weak pinning 
regime where the perturbation expansion is useful. 

Finally, the same methodology can be simply general- 
ized to a (d+ l)-dimensional geometry. The perturbation 
series diverges for d = 1 as amply demonstrated by vari- 
ous authors but useful results for all d > 1 can be 
obtained. In particular, results for the periodic potential 
H can be readily extended to the random case. 
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